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Abstract. We investigate the properties of motion in a 
map model derived from a galactic Hamiltonian made up 
of perturbed elliptic oscillators. The phase space portrait 
is obtained in all three different cases using the map and 
numerical integration of the equations of motion. Our nu- 
merical calculations suggest that the map describes very 
well the properties of regular motion in all three cases. 
There are cases where the map fails to describe chaotic 
motion and cases where the map describes satisfactorily 
the chaotic phase plane. We also compare the Lyapunov 
characteristic number and the spectra of orbits derived 
by the map and numerical integration in each case. The 
agreement of the outcomes is satisfactory. 
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1. Introduction 

The application of modern methods in the domain of 
Galactic Dynamics has been proved very fruitful during 
the last decade. Among them is the use of maps to describe 
galactic motion (Caranicolas 1990 , 1994 ), and the invari- 
ant spectra in galactic Hamiltonian systems (see Con- 
topoulos et al. |1995| , Patsis et al. p.997| ). 

Maps, derived from galactic type Hamiltonians, are 
useful in the study of galactic orbits because they are 
faster, in general, than numerical integration and allow 
a quick visualisation of the corresponding phase plane. 
Maps also can give the stability conditions for the peri- 
odic orbits. On the other hand, our experience based on 
previous work shows that the results given by the maps 
are in good agreement with those given by numerical inte- 
gration at least for small perturbations (see Caranicolas & 
Karanis 1999). On this basis, it seems that one has some 
good reasons to use a map for the study of galactic motion. 



In the present work we consider that the local (i.e. near 
an equilibrium point ) galactic motion is described by the 
potential 
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where u>i , u>2 are the unperturbed frequencies of oscillation 
along the x and y axis respectively, e > is the perturba- 
tion strength while a, /3 are parameters We shall study the 
case where u>\ — 102 =00. Without the loss of generality 
we can take u> = 1, that is the 1:1 resonance case. Then 
the Hamiltonian to the potential (1) is 
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where p x ,P v are the momenta per unit mass conjugate to 
x, y and h is the numerical value of H. 

Our aim is to study the various types of periodic or- 
bits, their stability and the kind of non periodic motion 
(regular or chaotic) in the Hamiltonian (2) for various val- 
ues of parameters a and (3 using the map corresponding 
to the Hamiltonian (2), as well as numerical integration. 
This resonance case is also known as the perturbed elliptic 
oscillators (Deprit 



1991 
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Innanen 1992, Caranicolas 1993). The x — p x Poincare 



phase plane derived by the map and the numerical inte- 
gration will be compared in each case. Of a special interest 
is the study of the chaotic motion. We shall try to find an 
answer to questions such as: 

1. Does the map describe in a satisfactory way the chaotic 
layers in the x — p x plane and, if so, how this behavior 
evolves by increasing e ? 

2. What are the differences, if any, in the Lyapunov Char- 
acteristic Number (LCN) found by the map and the 
numerical integration in the regular and the chaotic 
area? 

3. Are there any similarities in the invariant spectra de- 
rived by the map and the numerical integration? 
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The map and the stability conditions of the periodic 
orbits are given in Section 2. In the same Section we com- 
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Fig. 1. The x—p x (type A) phase plane derived by numer- 
ical integration. The values of the parameters are a = 1.2, 
P = 0.8, e = e esc = 0.03125, h = 2. 

pare the x — p x phase plane found by the map and numer- 
ical integration for some of the main different cases. In 
Section 3 we compare the LCNs and the spectra of orbits, 
derived using the map and numerical integration. Section 
4 is devoted to a discussion and the conclusions of this 
work. 



2. Map and stability conditions 

The averaged Hamiltonian corresponding to the Hamilto- 
nian (2) reads 



< H > = eaJ (A — J) [2 — cos 2</>] 



j 2 + (A - jy 



(3) 



Following the steps described in Caranicolas (1990) we 
find the map 



J n +i = Jn - 2eaJ„ + i(A - J n +x) sin2^ n 
4>n+i = ^„ + e(2a-3/3)(A-2J„ +1 ) 
—ea(A — 2 J n +i) cos 2</>„ 



(4) 



where A = h. The map describes the motion in the 
J — 4> plane and we return to x,p x variables through 
x = (2J) 1 / 2 cos(f>,p x = -(2J) 1 / 2 sin0. The fixed points 
of (4) are at 

(i) J — 0, for any tfi 

(ii) J = A/2, (j) = 0,tt (5) 
(in) J = A/2, cj) = ±tt/2 . 




Fig. 2. Same as Fig. |] derived by the map. 
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Fig. 3. The area near the center of Fig. [l] 



There are three distinguished case of the x — p x phase 
plane portrait in the Hamiltonian (2). Type A, when 
both fixed points (i) and (ii) are stable. Type B when 
points (i) are stable while fixed points (ii) are unstable. In 
type C phase plane, fixed points (i) and (ii) are unstable 
and stable respectively. Applying the stability conditions 



(see Lichtenberg & Lieberman 1983) we find after some 
straightforward calculations: Type A phase plane appears 
when a > [3 > a/3 (a > Q,/3 > 0). Type B appears when 
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Fig. 4. The x—p x (type B) phase plane derived by numer- 
ical integration. The values of the parameters are a = 0.2, 
(3 = 0.25, e = e esc = 0.125, h = 2. 

(3 > a (a > 0, (3 > 0). We observe the phase plane of type 
C when a > 3/3 (a > 0, /3 > 0) or if a < 0, /3 > 0. 

In all three cases there is, for a fixed value of the energy 
h, a value of the perturbation strenght e C sc, such as for 
e > e csc curves of zero velocity h— V = open and the test 
particle is free to escape. We do not consider cases where 
the curves of zero velocity are always closed that is cases 
where the Hamiltonian (1) has no e osc . The value of e GS c 
can be found usin g the method described in Caranicolas 
& Varvoglis ( 1984 ). For the type A phase plane we find 



l/[8h(a + 0)] , 



(6) 



while in the cases B and C e csc is given by the formula 



l/[16h/3] 



(J) 



Let us now go to see the three different types of the 
phase plane produced by the Hamiltonian system (2). In 
all numerical calculation we use the value h — 2. Fig. |] 
shows the type A x — p x phase plane derived by numeri- 
cal integration. The values of a, /3 are 1.2, 0.8 respectively 
while e = e csc = 0.03125. The motion is everywhere regu- 
lar except near the hyperbolic point in the center and in 
a thin strip along the separatrix.(see Fig. ||). Fig. ||is the 
corresponding figure produced by the map. As one can see 
the agreement is good. One significant difference is that 
the map is inadequate to produce the chaotic layer seen 
in Fig. |j. Also note that x max in Fig. [I] is greater than 2 
while in Fig. ^ is smaller. 

Fig. H and Fig. || show the type B phase plane derived 
using numerical integration and the map respectively. The 




Fig. 5. Same as FigM derived by the map. 



values of the parameters are a = 0.20, j3 — 0.25, while 
e = £csc = 0.125. The results are similar to those observed 
m Figs. and §. A gain the map describes well the real 
phase plane except in a small chaotic region near the two 
hyperbolic fixed points. Thus, our numerical calculations 
suggest that, the map (4) describes well the properties of 
motion in the Hamiltonian (2) up to the largest pertur- 
bation, that is e — e csc but it is insufficient to described 
the small chaotic region observed when using numerical 
integration. 

Type C phase plane is shown in Figs. || and 0. Fig. || 
comes from numerical integration while Fig. |^ was de- 
rived using the map. The values of the parameters are 
a = -1.20, /3 = 0.20, while e = e csc = 0.15625. This value 
of e was chosen to maximize the chaotic effects. The most 
important characteristic, observed in both Figs. |^ and ^ is 
the large unified chaotic region. In contrast with the two 
previous cases A and B, here the map reproduces satis- 
factorily the chaotic sea found by numerical integration. 
On the other hand, it is evident that the map describes 
qualitatively in a satisfactory way the areas of regular mo- 
tion around the elliptic points in the p x axis. Some differ- 
ences are observed in the area near the center between 
the two patterns. Another important characteristic, ob- 
served in case C, is that the chaotic areas are large even 
when e < e csc . Results, not shown here, suggest that con- 
siderable chaotic areas are observed in the phase plane 
derived using the map, when e > 0.08. Therefore we must 
admit that our numerical experiments show that, in the 
case when the Hamiltonian system (2) has small chaotic 
regions, the map is inadequate to produce them but it 
describes them satisfactorily when they are large. 
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Fig. 6. The x — p x (type C) phase plane derived by nu- 
merical integration. The values of the parameters are 
a = -1.2, f3 = 0.2, e = e esc = 0.15625, h = 2. 




Fig. 7. Same as FigJfl derived by the map. 



3. Dynamical spectra 

In this Section we shall study the spectra of the Hamilto- 
nian system (2) using numerical integration and the map 
(4). Before doing this, it is necessary to remember some 
useful notions and definitions. The "stretching number" 



(see Voglis & Contopoulos 1994, Contopoulos at al 



1995, Contopoulos & Voglis 1997) is defined as 



= In 
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where is the next image on the Poincare phase plane 
of an infinitesimal deviation & between two nearby orbits. 
The spectrum of the stretching numbers is their distribu- 
tion function 



S(ati 



AN(a) 
Nda 
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where A (a) is the number of stretching numbers in the 
interval (a, a + da) after N iterations. The maximal Lya- 
punov characteristic number can be written as 



LCN 



aS(a)da 



(10) 



in other words, the LCN is the average value of the stretch- 
ing number a. 

Today we know that the distribution of the successive 
stretching numbers forms a spectrum which is invariant 
with respect to (i) the initial conditions along an orbit 
and the direction of the deviation from this orbit and (ii) 
the initial conditions of the orbits belonging to the same 
chaotic region. 

In what follows we shall give the spectra S(a) of orbits 
of the Hamiltonian system (2) derived using the map (4) 
as well as numerical integration. In the case of the Hamil- 
tonian where t is a continuous quantity ( that is in the 
numerical integration of the equations of motion) we use 
for the derivation of the stretching numbers and spectra 
the method used by Contopoulos et al. (1995). All calcu- 
lations correspond to the case C because the results are 
much more interesting. 

Fig. H shows the spectra of two orbits (the first with 
solid line the other with dots) derived using numerical in- 
tegration. The orbits were started in the chaotic region 
with different initial conditions. It is evident that the two 
spectra are close to each other. Fig. || shows the spectra 
of the same two orbits found using the map. Again the 
two spectra are close to each other. The spectra in both 
cases were calculated for N — 10 6 periods. As one ob- 
serves the pattern shown in both Figs. || and || have the 
characteristics of the spectrum of orbits belonging to a 
chaotic region. Fig. 10 shows the LCNs for two chaotic 



orbits. Number 1 was derived using the map, while num- 
ber 2 was found using numerical integration. As one can 
see the mean exponential divergence of the two nearby 
chaotic trajectories, described by the map, is larger than 
that given by numerical integration. Nevertheless the two 
curves are qualitatively similar. 

In Figs, [ll] and 12 we give the spectra of the same reg- 
ular orbit derived by numerical integration and the map 
respectively. This is an orbit with initial conditions near 
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Fig. 8. The spectrum S(a) of two chaotic orbits derived 
using numerical integration. The values of the parameters 
are as in Fig. |^. Initial conditions x — 1.3, y — p x = and 
x = 1.1, y = 0, p x = 0.2. The value of p y > is found 
from the energy integral. 



the periodic orbit x = y = 0,p x = A 1 / 2 . The spectra were 
calculated for N = 10 6 periods. As one can see the agree- 
ment between the two spectra is good. From other cases 
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Fig. 10. LCNs for the same chaotic orbit given by the 
map 1 and numerical integration 2. Initial conditions x = 
1.3, y=p x = 0. 



we know that the agreement is much better if we calcu- 
late an orbit for N = 10 7 , or N — 10 8 periods. Both are 
"U-shaped", with two large and two small peaks. Such 
spectra are characteristic of quasi-periodic orbits, start- 
ing close to a stable periodic orbit (see Patsis et al. |1997|) . 
Fig. HH shows the LCNs for the orbit of Fig. [ll] derived 
using a map (dots) and numerical integration (solid line). 
Again one can see that the map describes well the quali- 
tative properties of regular orbits. 

Let us now come to the symmetry of the spectra. It 
was observed that the spectra of ordered orbits, starting 
close to a stable periodic orbit, are almost symmetric with 
respect to the a — axis, while, when we go far from the 
stable periodic orbit they become asymmetric. In order 
to give an estimate between closeness and symmetry we 
have made extensive numerical experiments near the exact 
periodic orbit x = 0, p x = h 1 / 2 = \/2. 

We define as "symmetry factor" the quantity 



1 N 



(11) 



As one can see, Q = corresponds to a perfectly symmet- 
ric spectrum, while Q = 1 to a totally asymmetric with 
all non zero values of S(a) to belong to either positive or 
negative a. The results Q vs p x for the map, in the case 
C, when e = 0.02, are shown in Fig. [l4|. Dots correspond 
to values given by the numerical experiments, while the 
solid line corresponds to the best fit 



Fig. 9. Same as Fig. U derived using the map. 
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Fig. 11. The spectrum S(a) of a regular orbit derived 
using numerical integration. The values of the parameters 
are as in Fig. [| Initial conditions x = y = 0, p x = 1.55. 
The value of p y > is found from the energy integral. 
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Fig. 12. Same as Fig. [Tl] derived using the map. 



where a = 0.2704, b = -0.7960 and c = 0.5986. Fig. || 
is the same as Fig. 14, for the Hamiltonian in the case 
C when e = 0.1. The best fit is now with a = 0.9100, 
b = —2.659 and c = 1.983. As we can see we have a 
second order polynomial growth of the asymmetry of the 
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Fig. 13. LCNs for the same regular orbit given by the 
map (dots) and numerical integration (solid line). Initial 
conditions as in Fig |ll|. 

spectrum as we move far from the stable periodic orbit. 
The value of N in all cases was 10 6 iterations. 



4. Discussion 

In this paper we have studied the regular and chaotic mo- 
tion in the Hamiltonian system (2) using a map, derived 
from the averaged Hamiltonian and numerical integration. 
Note that Hamiltonian (2) is integrable when a — (3 and 
when a — 3/3. Depending on the choice of the parame- 
ters a, p this Hamiltonian displays three different types 
of x — p x phase plane named as types A,B,C. In the case 
of the two first types of phase plane, the map describes 
in a very satisfactory way the real properties of orbits up 
to e = e csc . Extensive numerical calculations have shown 
that the map does not produce any chaotic regions in the 
cases A and B. Therefore, one concludes that the map fails 
to describe the small chaotic regions near the separatrix 
found by numerical integration. Here we must note that at 
the beginning we had started with smaller values of e, but 
we observed that the map persisted not to produce chaos 
up to e = e csc . This explains why e = e csc was chosen. 

The situation is quite different in the case of the phase 
plane of type C. Here the system has large chaotic regions, 
increasing when e approaches e C sc and the map describes 
them satisfactorily. Moreover, the comparison of the spec- 
tra and LCNs of orbits found using the map and numerical 
integration shows that one can trust the map. This is very 
important because the map is at least 10-20 times faster 
than numerical integration and one can make faster all 
the time consuming calculations. Note that in all cases 
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Fig. 14. Q vs p x for the map in case C, when e = 0.02 
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increases, with a second order polynomial dependence, as 
we are moving far from the stable periodic orbit. 

Finally the authors would like to make clear that (i) 
the proposed map has some qualitative similarities with 
the Poincare map of the original Hamiltonian system (2), 
although the numerical differences may be important and 
(ii) the results of this work correspond to the particular 
Hamiltonian system (2) and for the resonance case 1:1. 
Interesting results, on the spectra of orbits, in galactic 
Hamiltonians made up of harmonic oscillators in t he 4: 3 
resonance, have been given by Contopoulos et al. ( 1995 ). 
Also spectra of orbits have bee n stud ied, in the standard 
map, by Voglis & Contopoulos ( 1994 ), Contopoulos et al. 
( |1997| ). This paper was focused on the comparison of the 
results (structure of the phase plane and spectra of orbits) 
given by numerical integration and the map. 



Acknowledgements. The authors would like to thank the ref- 
eree Prof. G. Contopoulos for valuable suggestions and com- 
ments that helped to significantly improve the paper. 



References 

Caranicolas, N.D., 1990, Cel. Mech. 47, 87 
Caranicolas, N.D., 1993, A&A 267, 368 
Caranicolas, N.D., 1994, A&A 291, 754 
Caranicolas, N., Varvoglis, Ch., 1984, A&A 141, 383 
Caranicolas, N.D, Innanen, K.A., 1992, AJ 103, 1308 
Caranicolas, N.D. ,Karanis, G.I., 1999, A&A 342, 389 304, 374 
Contopoulos, G., Voglis, N., 1997,A&A 317, 73 
Contopoulos, C, Grousousakou, E., Voglis, N. 1995 A&A 
Contopoulos, C, Voglis, N., Efthymiopoulos, C, Froeschle, C, 

Gonczi, R., Lega, E., Dvorak, R., Lohinger, E., 1997, Cel. 

Mech. 67, 293 
Deprit, A., 1991, Cel. Mech. 51, 203 
Deprit, A., Elipe, A., 1991, Cel. Mech. 51, 227 
Lichtenberg, A. J., Lieberman, M. A., 1983, Regular and 

stochastic motion, Springer, Berlin Heidelberg, New York 
Patsis, P. A., Efthymiopoulos, C, Contopoulos, G., Voglis, N., 

1997, A&A 326, 493 
Voglis, N., Contopoulos, G., 1994, JPhys. A 27, 4899 



Fig. 15. Same as Fig. [14] for the Hamiltonian in case C 
when e = 0.1. 



the number of periods in the calculation of the spectra 
using the map or numerical integration was the same in 
order to be able to compare the corresponding results. 

We have also investigated the symmetry of the spec- 
tra of ordered orbits. It was found that the quantity Q 



